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Abstract 

The process of large-angle high energy electron-positron scattering with emission of 
one hard photon almost collinear to one of the charged particles momenta is considered. 
The differential cross section with radiative corrections due to emission of virtual and 
soft real photons calculated to a power accuracy is presented. Emission of two hard 
photons and total expressions for radiative correction are given in leading logarithmical 
approximation. The latter are illustrated by numeric estimates. A relation of the results 
with structure function formalism is discussed. 

PACS: 12.20.-m Quantum electrodynamics, 12.20.Ds Specific calculations [] 



1 Introduction 

The process of electron-positron scattering is commonly used for luminosity measurements at 
e + e~ colliders. It has almost pure electrodynamical nature and could therefore be described to 
any desired precision within a framework of perturbative QED. Nevertheless the accuracy of 
modern experiments is ahead of that provided by theory. A lot of work has recently been done 
to uplift the theoretical uncertainty to about one per mille under conditions of small-angle 
Bhabha scattering at LEP1 [|] and afterwards up to 0.05 — 0.06%@. 

The large-angle kinematics of Bhabha scattering process is extensively used for calibration 
purposes at e + e~ colliders of moderately high energies, such as 0, J /if), B, and c/r factories and 
LEP2. At the Born and one-loop levels the process was investigated in detail in || [|, |5], §], 01 > 
taking into account both QED and electroweak effects. 

In paper || we considered Bhabha scattering to O{o) order exactly improved by the 
structure function method. The latter, based on the renormalization group approach, allows 
to evaluate the leading radiative corrections to higher orders, including all the terms ~ (aL s ) n , 
n = 2,3, . . ., where L s = \n(s/m 2 ) is a large logarithm, s is the total center-of-mass (cms) 
energy of incoming particles squared and m is the mass of fermion. 
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To reach the one per mille accuracy it is required to take into account radiative correc- 
tions (RC) up to third order within the leading logarithmic approximation (LLA) and up to 
second order in the next-to-leading approximation (NLA). In a series of papers several sources 
of these corrections were considered in detail H [K], [11], [12]| . 

In a recent publication jTJ] the contribution due to virtual and soft photon corrections to 
large-angle radiative Bhabha scattering was calculated for the general case of hard photon 
emission at large angle with respect to all charged particles momenta. In the present work 
we are going to consider the complementary specific kinematics, in which the photon moves 
within a narrow cone of small opening angle 6 ^ 1 together with one of the incoming or 
outgoing charged particles. Thus, the result obtained here may be used in experiments with 
the tagging of scattered electron (positron) in detectors of small aperture 9 <C 1. 

Our paper is organized as follows. In Sec. 2 the Born level cross section of radiative Bhabha 
scattering is revised in the collinear kinematics of photon emission along initial (scattered) 
electron. We introduce here the physical gauge of real photon that is extensively used in 
the next sections. In Sec. 3 a set of crossing transformations which enables us to consider in 
some detail only the scattering type amplitudes of loop corrections to the process is described. 
Besides, we restrict ourselves to the kinematics of hard photon emission along initial electron. 
In Sec. 4 the corrections due to virtual and soft real photon emission in the case ki || p l are 
considered. The general expression for correction in the case of hard photon emission along 
scattered electron is given in Sec. 5. In Sec. 6 we consider a contribution (in LLA) coming 
from two hard photon emission and derive a general expression for radiative correction. In 
conclusion we discuss the relation with structure function approach and the accuracy of the 
results obtained. Some useful expressions for loop integrals are given in the Appendix and the 
results of numeric estimates are given in graphs. 



2 Born expressions in collinear kinematics 

Let us begin revising the radiative Bhabha scattering process 



e~(Pi) + e + (p 2 ) -> e-(pi) + e + {p' 2 ) + j(h) (1) 

at the tree level. We define the collinear kinematical domains as those in which the hard 
photon is emitted close (within a narrow cone with opening angle 9q <^ 1) to the incident 

(#i(2) = Pi(2)fei < #o) or the outgoing electron (positron) (9[, 2 ) — P'i(2)^i < #o) direction of 
motion. Because of the symmetry between electron and positron we may restrict ourselves to 
a consideration of only two collinear regions, which correspond to the emission of the photon 
along the electron momenta. The two remaining contributions to the differential cross section 
of the process (JT|) can be obtained by the substitution Q 



da. 



coll 



1 + Q 



Pl <-> P2 



P'2 



d(7 7 (fe 1 || ft) + d<7 7 (fci || p[ 



(2) 



To begin with, let us recall the known expression |13| in Born approximation for the general 
kinematics, i.e. assuming all the squares of the momenta transferred among fermions to be 
large compared to the electron mass squared: 
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s = (pi + p 2 ) 2 , si = (pi + pf,) 2 , t= (p 2 - p 2 ) 2 , h = (pi - pi) 2 , 
u = (pi - P 2 ) 2 > w i = (P2 - Pi) 2 , Xi = 2pifci, xi, 2 = Viol- 
in the collinear kinematical domain in which fci || p l the above formula takes the form 
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Here yj are the energy fractions of the scattered leptons and d<7 (pi(l — x),p 2 ) is the cross 
section of the elastic Bhabha scattering process. 

Throughout the paper we use the following relations among invariants 

Si + 1 + Ui = 4m 2 - xi ~ 0, s + ti + u = Am 2 + xi ~ 0. 



In the case k\ || pi we have 

<W(&i || pi) 
f 



^T— -(1 - a;)der (pi,P2), 



27T 2 Xi 

1 + (1 - x) 2 2m 2 



(6) 



x 



X'l 



These expressions could also be inferred by using the method of quasi-real electrons W$ and 
starting from the non-radiative Bhabha cross section. 
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After integration over a hard collinear (fei || p x ) photon angular phase space, the cross 
section of radiative Bhabha scattering in the Born approximation is found to be 
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where Lq = In (e9o/m) 



f3-3x + x 2 + 2cx(2 -x) + c 2 (l - x(l - x)Y 
[ (l-x)(l-c)a2 , 

And in the case k\ || p[ it reads 
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The simplest way to reproduce these results is to use the physical gauge for the real photon 
which in the beam cms sets the photon polarization vector to be a space-like 3-vector e\ having 
density matrix 



0. 



' fiv ^^T^U) \1> 
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These properties enable us to omit mass terms in the calculations of traces and, besides, 
to restrict ourselves to the consideration of singular terms (see Eq. (|10f)) only, both at the 
Born and one-loop level. As shown in |Tj|, this gauge is proved useful for a description of 
jet production in quantum chromodynamics; it is also very well suited to our case because 
it allows to simplify a lot the calculation with respect, for instance, to the Feynman gauge. 
What is more, it possesses another very attractive feature related with the structure of the 
correction to be mentioned below (see Appendix). 

With these tools at our disposal let us turn now to the main point. The contributions, 
which survive the limit 9q — > 0, arise from the terms containing 
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(10) 



Other omitted terms (in particular those which do not contain a factor Xi 1 ) can De safely 
neglected since they give a contribution of the order of 9$ which determines the accuracy of 
our calculations 



1 + 
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In the realistic case this corresponds to an accuracy of the order of per mille. 
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3 Crossing relations 



In this and the next section we shall consider the case ki \\ p 1 . In the case of photon emission 
along p[ one can get the desired expression by using the left-to-right permutation 





Figure 1: Some representatives of FD for radiative Bhabha scattering up to second order: (1 ) 
is the vertex insertion; (2) is the vacuum polarization insertion; graphs denoted by (3), (4) are 
of the L-type, (5) is of Gi-type, (6) is of G 2 -type, (7) is of B-type and (8) is of P-type. 



From now on we deal with scattering type amplitudes (FD) with the emission of hard 
photon by initial electron. This is possible due to the properties of the physical gauge. The 
contribution of annihilation type amplitudes may be derived by applying the momenta re- 
placement operation as follows: 

A|M|Lihilation = {QiPl <- -P2)} A|M| s 2 cat ten„g = {Ql} A I M I scattering • (13) 

In considering FD with two photons in the scattering channel (box FD) one may examine 
only those with uncrossed photons because a contribution of the others may be obtained by 
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the permutation p 2 <-» — p 2 . Thus the general answer becomes 



\M\] 



1 



»e{(l + Qi)[G + L] + — (1 + Q0(1 + g 2 )[ Sl t(5 + P)]}, 



with the permutation operators acting as 

QiF(s 1 ,t 1 , s, t) = F(t, s, t h st), Q 2 F(s, u, s lt u x ) = F(u, s, u h si). 



(14) 



4 Virtual and soft photon emission in k\ \\ p 1 kinematics 

One-loop QED RC (which are described by seventy two Feynman diagrams) can be classified 
out into the two gauge invariant subsets (see Fig.l): 

• single photon exchange between electron and positron lines (G,L-type); 

• double photon exchange between electron and positron lines (B,P-type). 

For L-type FD (see Fig. 1(3,4)) the initial spinor u(pi) is replaced by (a / (2tt)) A 2 k 1 eu(pi) , 
with 
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The relevant contribution to the matrix element squared and summed over spin states reads 
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The contribution of vertex insertion, vacuum polarization^] and Gi-type (see Fig. 1(1,2,5)) 
has the following form 
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Here A is as usual the IR cut-off parameter to be cancelled at the end of calculus against a 
soft photon contribution. 



For realistic applications one should also add to II the contributions due to fi and r leptons and hadrons. 



For the contribution of G 2 -type FD (see Fig. 1(6)) with four denominators we obtain 



A M U = 2 a 7T 
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(J - Jx)y + 2(2 z V (Jii - Ji + - xJ fe ) 
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It turns out that only the scalar integral and the coefficients before p±, k\ in the vector and 
tensor integrals give non-vanishing contribution in the limit 9 — > 



J 1 



d 4 k(l,k»,k»k u ) 
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in* (0)(1)(2)(9) 

(0) = k 2 - A 2 , (1) = k 2 - 2 Pl k, (2) = k 2 - 2p[k, (q) = k 2 - 2k( Pl - h) - X i, 
(aft)"" = a?b v + a"&", 

and the terms having no pi or £4 momentum in the decomposition have been omitted for their 
unimportance. 

The B-type FD shown in Fig. 1(7) with uncrossed legs gives 
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where the coefficients are associated with scalar, vector and tensor integrals over loop momen- 
tum 
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+ a, 22 \p 2 p' 2 y v . 

For P-type FD (see Fig. 1(8)) with uncrossed photon legs we have 
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Here we are using the definition (with tensor structures giving no contributions in the limit 
#o^O dropped) 



d A k (1, W) 
in 2 aoaia 2 a^a4 



(E, E lP 1 + E k kl EnP^ + E kk k»k\ + EtMK + pM), 



ao = k 2 — A 2 , a\ = k 2 — 2p\k, a 2 = k 2 — 2k(pi — k\) — xi, 
a 3 = k 2 + 2p 2 k, a 4 = k 2 — 2qk + t. 

Note that in the evaluating of P-type FD we are allowed to put k\ = xpi, thus keeping only 
Pi momentum containing terms in the decomposition. 
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Collecting all the contributions (for the explicit expressions of all the coefficients see 
Appendix) given above we arrive at the general expression for the virtual corrections with 
p = x[l + (e6/m) 2 ] <C s/m? 
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After integration over xi one gets additional large logs of the form L = L s + ln(^p/4). 
Terms containing the last factor have to be cancelled against a contribution coming from the 
emission of hard photon outside a narrow cone 9 < 8q 1 (and supplied by the same set of 
virtual and soft corrections), which was considered in fll~2| . In the case of two hard photon 
emission it is necessary to consider four kinematical regions, namely when both are emitted 
inside/outside a cone and one inside/another outside. 

Fortunately enough, the w-structure, which obviously violates factorization feature, does 
not contribute in LLA due to a cancellation of large logs in $. What for a correction to the 
above structure coming from P-type graph it vanishes in the sum of FD with crossed and 
uncrossed photon legs (for a more comprehensive account see Appendix). 

The total expression can be obtained by summing virtual photon emission corrections and 
those arising from the emission of additional soft photon with energy exceeding no Ae e. 

The emission of a soft photon is seen as a process factored out of a hard subprocess (in 
our case the latter is exactly a hard collinear photon emission) so this is seemingly come into 
an evident conflict with a hard collinear emission. Nevertheless, arguments similar to those 
given in the paper devoted to the problem of DIS with tagged photon |16| may be applied in 
the present paper: the factorization of the two in the differential cross section is present and 
we are, hence, allowed to consider a soft photon emission restricted as usual by 
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Thus the soft correction can be written as 
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where Mb denotes the matrix element of the hard photon emission at the Born level and in 
the kinematics ki \\ p 1 it reads 
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Now let us check the cancellation of the terms containing L\. Indeed it takes place in the sum 
of contributions arising from emission of virtual and soft real photons. To show that we bring 
the soft correction into the form 
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where q are the cosines of emission angles of i-th particle with respect to the beam direction 
(pi in cms), cy 2 > is the cosine of the angle between scattered fermions in cms of the colliding 
particles and yi are their energy fractions and in the case fci || p l we have 



1 2 yi y 2 ' 2 2y 2 {l-x) 1 ; 

Then the cancellation of infrared singularities in the sum is evident from comparison of 
Eqs. (|20|j2^) . The terms with ln(Ae/e) should be cancelled when adding a contribution of a 
second hard photon having energy above the registration threshold As. 

The complete expression for the correction in the case ki \\ p 1 reads 
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5 Kinematics k\ \\ p[ 



We put here a set of replacements one can use in order to obtain the modulus of matrix 
element squared and summed over spin states for the case ki || p[, starting from the analogous 
expression for ki \\ p 1 (Eq. (|T4])) and using the replacement of momenta p\ <-> —p[,P2 ^ — P^- 
The last operation results in the following substitutions: 
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Then under these permutations the expression for virtual corrections given in Eq. 
transformed yielding the following result for the collinear domain k\ || p[ 
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and $, $ derived upon applying a set of replacements from Eq. (P7D on the quantities <3>, 
The contribution from the soft photon emission is described by 
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The total correction for the case ki || p[ has the following form 
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Performing the integration over a hard photon angular phase space (inside a narrow cones) 
we put the RC to the cross section coming from virtual and soft real additional photons valid 
to a logarithmic accuracy in the form 
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In the Fig. 2(a,b) given are the ratio of S/(L t Si) versus x for the collinear kinematics consid- 
ered above. 



6 Two hard photon emission and results in LLA 



Turning to the structure of the result obtained, it should be noted that all the terms quadratic 
in large logarithms L tl ~ L S1 ~ L u ^> L p are mutually cancelled out as it should be. 

From the formula ( p6|) it immediately follows that (upon doing an integration over a hard 
photon angular (within a narrow cone) phase space) the w-term that is not proportional to T, 
which is in fact the kernel of the non-singlet electron structure function, is not dangerous in 
the sense of a feasible violation of the expected Drell-Yan form of the cross section, because 
it does contribute only at next-to-leading order. 

Performing the above mentioned integration and confining ourselves to LLA we get for the 
sum of virtual and soft photons 



d(T 7(S+V) 

dxdc 



do"o ol 
dxdc Ti 



4 In 1 ln(l 

e 3 2 v 



x) - \n{y x y 2 ) 



(32) 



The LLA contribution coming from the emission of second hard photon with total energy 
exceeding Ae consists of a part corresponding to the case in which both hard photons (with 
total energy ex) are emitted by initial electron 



dxdc 



d<jQ a 
/ . 

dxdc 7T 

■ 1 + z 2 
.1-z 



xVg\l 



4(1 + (1 -x) 
21n(l - z) 



+ - ln(l - x) - ln — 



3 

mz H — 
2 



l + z 
H In z 



3 
' 4 

l + z 



(33) 



and the remaining part which describes the emission of second hard photon along scattered 
electron and positrons. The latter, upon combining with the part of contributions of soft and 
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virtual photons to our process 



dcx 7 3 a 



-L 



' Ae 3 
In — + - 

e 4 



dxdc 7T 

may be represented via electron structure function in the spirit of the Drell-Yan approach 



/ d(J v 

Mxdc ; 



a 1 



x 



27T X 
do- (pi(l -x),^2P2jgl,g2 

dc 



L / dz 2 dz 3 dz 4 T>(z 2 )V(z 3 )V(z 4: ) 



(34) 



with the non-singlet structure function JHJ 



Viz) 



5(1 



2tt 



a 



2tt 



2! 



(35) 



p (1 ' 2) 0?) = hm - ^)pi 1,2) + e(i - A-^p^^j, 



21nA 



1 + z 2 



p 



(2) 



2 In A + 



2tt 2 



2 ' ° v ' 1-2 

These functions describe the emission of (real and virtual) photons both by final electron 
and by positrons. The multiplier before the integral stands for the emission of an hard photon 
by the initial electron. Thus Eq. ([$4]) actually represents the partially integrated Drell-Yan 
form of the cross section. Quite the same arguments are applicable to the second case in which 
an hard photon is emitted by the final electron. 

The cross section of the hard sub-process e{p\Z\) + e(p 2 z 2 ) — > e(qi) +e(q 2 ) entering Eq. (|34l ) 
has the form 

' • (36) 



da (zipi,z 2 p 2 ;qi,q 2 ) 87ra 2 



dc 



Z l + Z \ + Z 1 Z 2 + 2c(z\ - z\) + (?{z\ +Z 2 - ZiZ 2 ) 



Z X (\ - C)(«i + Z 2 + C(^ 2 - Zx)) 2 

The momenta of scattered electron qi and positron q 2 are completely determined by the 
energy-momentum conservation law 

2z 1 z 2 



£- 



Z\+ z 2 + c(z 2 - Z\) 
c = cos <hTp 1 , zi sin q^p 1 



q°i+q 2 =e(z 1 + z 2 ), 

= z 2 wxq^p^. 



In general their energies differ from those detected in experiment e[,e 2 , namely 



Ql z 3, 



q 2 Zi, 



whereas the emission angles are the same in LLA. 



Collecting the two expressions presented in Eqs. (p2| , |33|) one can rewrite the result in LLA 



as 



da^ 



dxdc 



Si 



da2 

dxdc 



fcillPi 
-1 + 



( da 7 \ 
\dxdc) killpi 

2 



{l + 5i}. 



L 



fcillPi 



3 Kyyj 4 1+ 1 



x 



(37) 
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For the case ki || p[ the correction is found to be 



da^< 



dxdc 



/ d ^o v 
Mxdc 



fexllPi 

5v 

fcillPi 



(dxdc) fci||p , {1 + 5l ' } 



/ d<xS v 



1 + -L 

7T 



X 



dxdc / fciiip; 
a l + (l-x) 2 

27T X 

do- (zipi,z 2 p2;gi,g2) 
dc 



2 + x^ 2) (l 



3 4(l + (l-x) 2 ) 



L' J dz 1 dz 2 dz i V(z 1 )V(z 2 )T>(z A ) 



(38) 



3S 



Pi 



with L' = L + 2 ln(l - x) . 

For the case when the energies of scattered fermions are not detected the expressions ([J 
may be simplified due to / dzV(z) = 1 and z^, ^-independence of the integrand in fci 
kinematics (^-independence in ki || p[ case). 

The x-dependence of 8% are shown in the Fig. 3 for different values of the cosine of scattering 
angle c. For a hard photon emission by final particles the correction S[ strongly depends on the 
experimental conditions of particles detection: the energy thresholds of detection of scattered 
fermions. This dependence for S\ is much more weaker, namely about 1%. 

In conclusion let us recapitulate the results given in Eqs. ( |37l , |3"8"l) . They both respect the 
Drell-Yan form for a cross section in LLA. Nevertheless a certain deviation away from RG 
structure function representation at a second order of PT in k\ || p 1 kinematics is observed. 
The term destroying expectations based on RG approach comes from definite contribution of 
a soft photon emission, the term with \n(yiy 2 ) in Eq. ( |3?| ) which cannot be included into the 
structure function approach. Its appearance is presumably a mere consequence of a complicate 
kinematics of 2 — > 3 type hard subprocess (see WW)] for such a kind of processes the validity 
of the Drell-Yan form for a cross section was not proved so far. Another possible way out is 
a careful analysis of a conflict between a soft and hard collinear photon emission. We have 
used the factorized form of a soft photon emission (|22|) under the condition (|21~D. But, to the 
moment, this representation in the peculiar case at hand is not rigorously proved as well. 

The accuracy of our calculations of virtual and soft photon corrections is determined by 
the omitted terms of the order of 



V 7T S 7T 



(39) 



which corresponds to a per mille level. The accuracy of the correction coming from two hard 
photon emission is determined by 0((a/ir) ln(4/^)) and at 1% level. 
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Appendix 

Here we give the expressions for the quantities associated with G-type integrals: 



J 



■h 



Jk 



hi 



J 



lk 



1 



Xiti 



7T~ 



-2L x L tl + 2L tl L p - L t - 2Li 2 (x) - — 



dz In z 



hXi J l-zl-Xz hxi 

n 

dz z\nz 



A ( x 
1 + 



p — x 



hXiPJ l-zl-Xz 
dz 



tiXiJ (l-z)(l-Az) 

n 

zdz 



z In z \ 
1 + - I , 



1 + 



l-z / 
zliaz\ 



A = Li 2 (l-p) 



TV 



6 

A + $ 
hXi 



Li 2 (x) + L p ln(l — x), X 



x 



P 



Xi 

P = — 2 ' 



(A.1) 



In the limit p > 1 we have 



$ = xiM + hxi(Ju ~h + xJik - xJ k ) = -- + 0{p x ) 

and that is the reason why w-structure does contribute only to next-to-leading terms. 

In general the expression for 5-denominator one-loop scalar, vector and tensor integrals 
are some complicate functions of five independent kinematical invariants (in the derivation we 
extensively use the technique developed in |Tj|). In the limit m 2 <C xi s ~ — t they may 
be considerably simplified because of singular 1/xi terms only kept: 



E 
Ex 

-D0124 

KeL> 234 



— -D0124 + 7-D0123? 

Si t 



-xE k 
1 

1 



2X3 



(-D0134 - (1 - £)A,234 - xD 123 4 + XlE) 



Li + 2LJn 



X 



1 — X 



In 2 



x 



2tt 2 



SX\ 
1 



71 



1 — X 
2 



Lg - 2L S1 L P - 2L S L\ + — + 2Li 2 (x) 

o 



-^si + 2L Sl L\ — 2L p L sl + 2L sl L t 



5tt 2 



(A.2) 
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Ke-D 134 
3?eDi234 



1 

si 



L\ + 2L S L\ - 2(L tl + \n(x))L s + 2L s L t 
1 -Lg + 2L s (L tl + ln(x)) + 2L Sl L> 



S\Xt\ 



111 
6 

7tt 2 



6 



The structure En + xE\k has the form l/(sxi)/(^, Xi) an d wm vanish after performing the 



operation (l + Q 2 )sitP given in (TA) which yields a contribution of P-type graphs with crossed 
and uncrossed photon legs. 

The following coefficient for the scalar integral is obtained in the calculation of B-type FD: 



B 



1 

Sit 



^si + %L Sl L\ — 2L Sl L p + 2L Sl L t + 



7T 



(A.3) 



For the vector integral coefficients we get 
1 



-n 2 s 1 + 2«iLi 2 (l - p) - SiL; + tL 2 Si - 2tL Sl L t 



1 



2s x t 
1 



2tt 2 



2Li 2 (l - p) - 2L\ + AL S1 L P - 2L ai L, 



(A.4) 



7T 



2«iLi 2 (l - p) + — (4«i + 6t) + (t - 2«i)L* - SiL t " + A Ul L Sl L p + 2 Sl L Sl L t 
o 



2si«it 

The relevant quantities for tensor B-type integrals are 



di>2' 

a V2 
Jo 



P 



S\t \p — 1 

1 



Lp-L t j, a g = -^-[(L S1 - L t ) 2 + n 2 ], 



P 



, 2 [(L t -L Sl f + n 2 } + ^-(L sl -L t )- — 
2u{ tu\ sit \ p — 1 



1 

si 



3 4 2 

-L^-2L Sl L p -Li 2 (l-p)-^- 



(A.5) 



As has been mentioned in the text, the physical gauge exploited provides a direct extraction 
of the kernel of the structure function out of the traces both in the tree- and loop-level 
amplitudes. The pattern emerging 



(Pi - h + m)e(pi + m)e(pi - k\ + m) = 4(pie) 2 (pi - ki) - (?X\k\ 

7H (1 — x)Ypi, 

kiHpi + rn)e(pi — ki + m) « (1 — x) (2 -W — y\ p\ 

V 1 - x J 



(A.6) 



shows this clearly. 
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Figure captions 



Figure 2. 

The ratio ^ (see Eq. ([H])) versus x = — for the case: 

a) fei || p x . 

b) fci || p[. 



Figure 3. 

The x-dependence of 5\ (see Eq. fl57p) 

for different values of the cosine of scattering angle c. 



Other parameters chosen are: 9q = 0.1, £ = 1 GeV. 
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